Renormalization group study of the four-body problem 
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We perform a renormalization group analysis of the non-relativistic four-boson problem by means 
of a simple model with pointlike three- and four-body interactions. We investigate in particular the 
region where the scattering length is infinite and all energies are close to the atom threshold. We 
find that the four-body problem behaves truly universally, independent of any four-body parameter. 
Our findings confirm the recent conjectures of Platter et al. and von Stecher et al. [H-Q that the 
four-body problem is universal, now also from a renormalization group perspective. We calculate 
the corresponding relations between the four- and three-body bound states, as well as the full bound 
state spectrum and comment on the influence of effective range corrections. 
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I. INTRODUCTION 

During the last decade few-body physics experienced 
a renewed interest due to the advent of experiments with 
ultracold atomic gases. Whereas the study of few-body 
physics in nuclear systems is hindered by the large com- 
plexity of the interparticle potentials, the interactions in 
ultracold atomic gases are describable to high accuracy 
with very simple short-range models. In addition, ul- 
tracold atomic gases become even more attractive as an 
ideal theoretical and experimental playground since they 
do not only offer excellent experimental control but also 
the amazing possibility of tuning the two-body interac- 
tion strength over a wide range using so-called Fcshbach 
resonances 

This made it possible, that.about forty years after V. 
Efimov's seminal prediction [5[ of the existence of univer- 
sal three-body bound states in systems with large two- 
body interactions, first evidence in favor of the presence 
of these states had been found in the remarkable experi- 
ment by Kraemer ct al. in 2006 Q. In his work, Efimov 
predicted the existence of infinitely many trimer states 
for infinitely large scattering length where the two-body 
interaction is just on the verge of having a bound state. 
The energy levels form a geometric spectrum and the 
three-body system is found to be universal in the sense 
that apart from the s-wavc scattering length a only one 
piece of information about the three-body system enters 
in the form of a so-called three-body parameter Q. The 
findings of Kraemer et al. stimulated extensive activity 
in the field of three-body physics, both experimentally 
[c!Jl2T| and theoretically; for recent reviews on also the 
latter see @, EH and references therein. As a result, the 
Efimov effect in three-body systems is a well-understood 
phenomenon today. 



The next natural step is to raise the question what the 
physics of four interacting particles may be. Early at- 
tempts towards an understanding of this system were 
made in the context of nuclear physics using a variety 
of approaches [T3 - [l7j . Also the four-body physics of 
4 He atoms has been investigated in much detail, for an 
overview see, e. g. [Tc| . The simpler four-body physics 
of fermions with two spin states, relevant for the dimcr- 
dimer repulsion, has also been studied fl9| . 
In their pioneering work, Platter and Hammer, et al. 
[3, [|| investigated the four-body problem using effective 
interaction potentials and made the conjecture that the 
four-boson system exhibits universal behavior. They also 
found that no four-body parameter is needed for a self- 
consistent renormalization of the theory. Calculating the 
energy spectrum of the lowest bound states in depen- 
dence on the scattering length a the existence of two 
tetramer (four-body bound) states associated with each 
trimer was conjectured. 

Recently, von Stecher, D'Incao, and Greene [H, [2(| in- 
vestigated the four-body problem in a remarkable quan- 
tum mechanical calculation. They found that the Efimov 
trimer and tetramer states always appear as sets of states 
with two tetramers associated with each of the trimer lev- 
els and calculated the bound state energy spectrum of the 
lowest few sets of states. The calculation suggests that 
the energy levels within one set of states are related to 
each other by universal ratios, which were obtained from 
the behavior of these lowest sets of states. In accordance 
with the results of Platter et al. [H, 0] the absence of 
any four-body parameter was also demonstrated. In or- 
der to find experimental evidence of the tetramer states 
extremely precise measurements are required. Remark- 
ably, Ferlaino et al. were able to observe signatures of the 
lowest two of the tetramer states in a recent experiment 
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While the calculations by Platter ct al. [l], Q and von 
Stecher et al. 0, [13] rely on quantum mechanical ap- 
proaches, in this work we want to shed light onto the 
four-body problem from a different perspective. A lot of 
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FIG. 1: (Color online) The generalized Efimov plot for four 
identical bosons. Here, we plot the energy levels of the vari- 
ous bound states as a function of the inverse s-wave scattering 
length a as numerically calculated in our approximative, ef- 
fective theory. In order to improve the visibility of the energy 
levels we rescale both the dimensionless energy E/A 2 and the 
dimensionless inverse scattering length a~ 1 /A where A de- 
notes the UV cutoff of our model. Also, we only show the 
first three sets of Efimov levels. The solid black line denotes 
the atom-atom-dimer threshold, while the dotted black line 
gives the dimer-dimer threshold. In the three-body sector 
one finds the well known spectrum of infinitely many Efimov 
trimer states (green, dashed) which accumulate at the unitar- 
ity point E^ = a -1 = 0, indicated by the orange star. In our 
pointlike approximation the four-body sector features a single 
tetramer (solid, red) associated with each trimer state. 



insight into the three-body problem had been gained from 
effective field theoryand renormalization group (RG) 
methods @, EH, [H, and it is desirable to apply these 
also to the four-body problem. In this paper we will make 
a first step towards such a description that is complemen- 
tary to the previous quantum mechanical approaches. 

Of special interest is the further investigation of uni- 
versality in the four-body system. In this context the 
so-called unitarity point, illustrated by the star in Fig. 
[TJ is of particular importance. In this limit not only the 
scattering length a is infinite but also all binding energies 
in the problem accumulate at the atom threshold at zero 
energy. Only at the unitarity point physics becomes truly 
universal in the sense that for example the ratio between 
the binding energies of consecutive trimer levels assumes 
exactly its universal value, E n+ i/E n = exp(— 2-k/sq), 
with so ~ 1.00624 the so-called Efimov parameter. The 
unitarity point is therefore the most interesting one from 
a theory point of view. Unfortunately, in the previous 
calculations only few lowest lying states were determined. 
The major advantage of the present RG approach is that 
it allows to investigate analytically the complete spec- 
trum and to address directly the unitarity point in order 
to extract the universal relations between the three- and 
four-body bound states in this limit. 
We investigate the renormalization group behavior of the 
relevant four-body interactions with an approximate, but 



simple and physically intuitive model which allows only 
for pointlike (i.e. momentum independent 1 ) three- and 
four-body interactions. In the three-body problem uni- 
versality manifests itself in an RG limit cycle of the three- 
body coupling. We find that this three-body limit cycle 
leads in turn to a "self-sustained" limit cycle of the four- 
body sector leaving no room for any four-body parame- 
ter. 

The RG method allows furthermore for computations 
away from the unitarity point. We calculate the bound 
state energy spectrum in the pointlike approximation 
(see Fig. [T|) and investigate how the relations between 
tetramer and trimer states approach the universal limit 
as one comes closer to the unitarity point. 
The paper is structured as follows. In Sec. [XT] wc 
introduce the functional renormalization group (FRG) 
method and set up the microscopic model. Sees. IIIII 
and IIVI are devoted to the FRG analysis of the two- and 
three-body sector. In Sec. [V] we discuss the four-body 
sector and present our numerical results. Our findings 
are summarized in Sec. IVII 



II. METHOD AND DEFINITION OF THE 
MODEL 

In this work we are interested in the computation of 
the few-body properties, such as the bound state spec- 
trum, of four identical bosons. In a quantum field theory 
approach the information about these properties can be 
extracted from the effective action T which is the generat- 
ing functional of one-particle irreducible vertex functions 
r'"' and which contains all information about a given 
system. The computation of T is a very complicated task 
as quantum (and, as in the case of nonzero density, statis- 
tical) fluctuations have to be integrated out on all length 
and therefore momentum scales q. In order to cope with 
this task we rely on the functional renormalization group 
[24| , for detailed reviews we refer to [III [26[ . 
The central quantity of the FRG is a scale dependent 
effective action functional, the so-called effective flowing 
action T^. The effective flowing action T^, which in- 
cludes all fluctuations with momenta q > k, interpolates 
between the classical action S at some ultraviolet (UV) 
cutoff scale k = A and the full quantum effective action 
r in the limit k — > 0. The underlying idea is similar 
to Wilson's idea of momentum shell-wise integration of 
fluctuations. The evolution of Tk is governed by the Wet- 
terich equation [24| . which is an exact, non-perturbative 
RG equation. It reads 

d k T k = ^Tr (if } + i? fc ) ~* d k R k , (1) 



Throughout the paper we will use the term momentum inde- 
pendence to refer to combined spatial momentum and frequency 
independence. 
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where is the flowing, full inverse propagator and 
the trace Tr sums over momentum q and Matsubara fre- 
quency q as well as the internal degrees of freedom such 
as species of fields. The dependence on the RG scale k is 
introduced by the regulator Rk- At the UV scale k = A 
the effective flowing action Tk equals the classical action 
S and as we want to consider dilute atomic gases the UV 
scale A is set to be of the order of the inverse Bohr radius 
Oq . For most problems, quantum and statistical fluctu- 
ations will generate infinitely many terms in Tk- Due 
to this fact it is in practice impossible to solve Eq. ^ 
exactly. Therefore one has to decide for a truncation of 
Tfc , which in turn corresponds to solving the theory only 
approximately. 

In this work we investigate the four-boson problem by 
approximately solving Eq. ([I]). Our truncation for the 
Euclidean flowing action is given by a simple two-channel 
model 

r k = f{tp*(d T -A + E^)tp 

J x 

+ /3((j)*(j)*(jyipip + <txjxt>* ip*ip*) +jcj)*tp*ip*(tyiptp}, (2) 

where A denotes the Laplace operator and we use the 
natural, non-relativistic convention 2M = h = 1 with 
the atom mass M. ip denotes the field of the elementary 
bosonic atom, while the dimer, the bosonic bound state 
consisting of two elementary atoms, is represented by the 
field </> ~ Both the atom and the dimer field are sup- 
plemented with non-relativistic propagators with energy 
gaps Eip and m^, respectively. In our approximation the 
fundamental four-boson interaction ~ X^(ip*ip) 2 is medi- 
ated by a dimer exchange, which yields A^, = — /i 2 /m| in 
the limit of pointlike two-body interactions. The dynam- 
ical dimer field cf> allows us to capture essential details 
of the momentum dependence of the two-body interac- 
tion. We introduce a wave function renormalization fac- 
tor Aif, for the dimer field in order to take into account 
an anomalous dimension of the dimer field <j). The only 
nonzero interaction, present at the microscopic UV scale 
k = A, is taken to be the Yukawa-type term with the cou- 
pling h. Together with the microscopic value of the 
Yukawa interaction h at the UV scale can be connected 
to the effective range r c ff in an effective range expansion. 
The atom-dimer interaction Xad as w °h as the various 
four-body interactions A^, /3, and 7 vanish at the UV 
scale and are built up via quantum fluctuations during 
the RG flow. 

At this stage we want to emphasize the meaning of the 
term pointlike approximation which must not be con- 
fused with the notion of a zero-range (contact) model. 
Consider for example the two-body contact interaction 
~ \,p(ip*ip) 2 , which has no momentum dependence on 
the UV scale, k = A. In order to describe the scattering 
of two particles in quantum mechanics one proceeds by 



solving the two-body Schrodingcr equation. From this 
one obtains the well-known result for the zero-range s- 
wave scattering amplitude fo(p) = (—a 1 — ip)~ l (with 
p = \p\ denoting the momentum of the colliding par- 
ticles). The scattering amplitude becomes momentum 
dependent. In the RG approach one deals with the effec- 
tive vertex A^, which varies with the RG scale k. On the 
UV (k = A) scale is momentum independent. When 
including more and more quantum fluctuations - mean- 
ing lowering the RG scale k from A to eventually k = 
- the effective vertex function A^, assumes a momentum 
dependence which in the IR limit k = is equivalent to 
the result for /o in the zero-range model. In a pointlike 
approximation one ignores this generated momentum de- 
pendence. In the simple model Eq. the three- and 
four-body sector is treated strictly in the pointlike ap- 
proximation. However, in the two-body sector the mo- 
mentum dependence of effective vertex A^, is captured 
by the exchange of the dynamic (i.e. momentum depen- 
dent) dimer propagator, such that the two-body sector is 
treated beyond the pointlike approximation. 
In the general case of nonzero density and temperature 
one works in the Matsubara formalism and the integral in 
Eq. ([2]) sums over homogenous three-dimensional space 

and over imaginary time J x = J d 3 x dr. Although 
our method allows us to tackle a full, many-body prob- 
lem at finite temperature in this way, we are interested 
solely in the few-body (vacuum) physics in this paper, for 
which density n and temperature T vanish. For T — 0, J 
reduces to an integral over infinite space and time. Our 
truncation (J5J) is based on the simple structure of the non- 
relativistic vacuum and, as demonstrated in (23l.[27|. nu- 
merous simplifications occur when solving Eq. (JTJ com- 
pared with the general, many-body case. The flowing ac- 
tion {2J has a global U(l) symmetry which corresponds 
to particle number conservation. In the vacuum limit it 
is also invariant under spacetime Galilei transformations 
which restricts the form of the non-relativistic propaga- 
tors to be functions of d T — A for the atoms and d T — A/2 
for the dimers. All couplings present in Eq. ^ are al- 
lowed to flow during the RG evolution and are taken to 
be momentum-independent in Fourier space as explained 
above. 

Besides the ansatz of Tk we must choose a suitable reg- 
ulator function Rk in order to solve Eq. ([T]). Based on 
our recent treatment of the closely related thrce-fermion 
problem [28|, [29[, we choose optimized regulators, 

= (k 2 -q 2 )e(k 2 -q 2 ), 
R* - ^{k 2 -q 2 )e{k 2 -q 2 ), (3) 

wi th q = \q\. These regulators are optimized in the sense 
of [Ifl [3(| and allow to obtain analytical results. 
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III. TWO-BODY SECTOR 

A remarkable and very useful feature of the vacuum 
flow equations is comprised by a special hierarchy: the 
flow equations of the iV-body sector do not influence the 
renormalization group flows of the lower N — 1-body sec- 
tor [23| . For this reason the different iV-body sectors can 
be solved subsequently In this spirit we first solve the 
two-body sector, then investigate the three-body sector 
in order to finally approach the four-body problem. 

The solution for the two-body sector can be found ana- 
lytically in our approximation 2 (for the analogous prob- 
lem considering fermions, see [28l [3lT|). The only run- 
ning couplings in the two-body sector are the dimer gap 
and its wave function renormalization A$. The flow 
equations of the two-body sector are shown in terms of 
Feynman diagrams in Fig. [2ja) and read 



dtA* 



h 2 



k 5 



12^ 2 (fc 2 + ^) 2 ' 
h 2 k 5 
12tt 2 (£; 2 +£^) 3 ' 



(4) 



where t = In 4 . As there are no possible nonzero flow 
diagrams for the Yukawa coupling h, it does not flow in 
the vacuum limit. 

The infrared (IR) values of the couplings h and can 
be related to the low-energy s-wave scattering length a 
via 



the closed channel with respect to the atom-atom thresh- 
old HH and the coupling h is proportional to the width of 
the associated Feshbach resonance being a function of the 
strength of the coupling to the closed channel. The choice 
A<p(A) = 1 then corresponds to the so-called characteris- 
tic length r* = — |r e ff often used in literature fH l33j. 
In the limit of large, positive scattering length there 
exists a universal, weakly bound dimer state. In or- 
der to find its binding energy we calculate the pole of 
the dimer propagator, corresponding to the condition 
m\(E$, k = 0) = 0, which yields in the limit £^,/A 2 « 1 



Er 



-2E, b = -2 




1 - 



2r, 



off 



(7) 



In the limit h — > 00, corresponding to r e s — > one re- 
covers the well-known result Ed = —2/a 2 . The dimer 
bound state energy is shown as a function of the inverse 
scattering length in Fig. Q] (black solid line). The de- 
viation from the universal 1/a 2 scaling for large inverse 
scattering lengths is due to the finite size of h which is 
taken to be h 2 /A = 10 in Fig. [TJ In the regime of small 
scattering length a one finds a crossover of the behavior 
of the dimer binding energy which then has the limiting 
behavior Ed = 4/(ar e ff). 



h 2 (k = 0) 



167rm 2 (fc = 0,£</> = 0) 



(5) IV. THREE-BODY SECTOR 



Knowing the analytical solution of the two-body sec- 
tor, this relation can be used to fix the initial values of 



our model, 
find 



For the UV value of the dimer gap m , we 



h 2 

(A) = -— a 
16tt 



12tt 



The first term fixes the s-wave scattering length ac- 
cording to Eq. ([5]), while the second term represents a 
counterterm taking care of the UV renormalization of 
the two-body sector. Finally, the last term accounts for 
the fact that the dimer consists of two elementary atoms. 
Additionally, we choose A<j,(A) = 1 which corresponds 
to the effective range r c g 



647T 



The action @ can also be used for a quite accurate 
description of Feshbach resonances. In this context Eq. 
@ is referred to as resonance model. In such a model m 2 ^ 
is proportional to the detuning energy of the molecule in 



2 Remarkably, the atom inverse propagator (one-body sector) is 
not rcnormalized in the non-rclativistic vacuum. 



The bound state spectrum of the three-body sector is 
much richer than the one of the two-body system. In 
his seminal papers [f| Efimov showed the existence of 
an infinite series of three-body bound states for strong 
two-body interactions. These energy levels exhibit a 
universal geometric scaling law as one approaches the 



(6) unitarity point £y 



0. Remarkably, these 



three-body bound trimer states exist even for negative 
scattering lengths a where no two-body bound state is 
present; they become degenerate with the three-atom 
threshold for negative scattering length and merge into 
the atom-dimer threshold for positive a. An additional 
three-body parameter is needed in order to determine 
the actual positions of the degeneracies [22[ . In this 
section we want to shortly review how Efimov physics 
can be treated within our approach. For a more detailed 
account on that matter and an application to the 
three-component 6 Li Fermi gas we refer to [23l. [28l. [29j . 

In our truncation, the three-body sector contains a sin- 
gle, pointlike (j)*ip*(j)ip term with a coupling Xad, which 
is assumed to vanish in the UV. It is build up by quan- 
tum fluctuations during the RG flow and the correspond- 
ing Feynman diagrams of the flow equation for Xad is 
shown in Fig. [U(b). First we investigate the unitarity 
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point, E^, = a -1 = 0. For this limit we are able to ob- 
tain an analytical solution for the flow equation of Xad 
while away from unitarity we have to rely on a numerical 
solution. 

At the unitarity point all intrinsic length scales drop 
out of the problem and the system becomes classically 
scale invariant. At unitarity, the Yukawa coupling h is di- 
mensionlcss and the only (extrinsic) length scale present 
is the inverse ultraviolet cutoff A -1 , which defines the 
validity limit of our effective theory. 
In our approximation the dimer field (f> develops a large 
anomalous dimension r\ = — ®±A± = 1 at unitarity which 
is consistent with the exact solution of the two-body sec- 
tor [23|, [I?]] ■ As the atom and dimer propagators have 
vanishing gaps in the IR the two-body sector respects a 
continuous scaling symmetry. 

In order to find the solution of the three-body sec- 
tor we switch to the rescaled, dimensionless coupling 
Xad = j^Xad- One finds that the flow equation for 
Xad becomes independent of k and h, 

+ itM ■ 

c 

As was demonstrated in [23|, [HI , the behavior of the so- 
lution of this type of flow equation is determined by the 
sign of the discriminant D of the right hand side of Eq. 
dSJ which is D = b 2 - 4ac < 0. Eq. (|5J) can be solved 
analytically and one finds 

-b + x/^Dtan (^-(t + 5)) 
Xad (t) = ^ L , (9) 

where S is connected to the three-body parameter and 
determines the initial condition. Most remarkably, the 
three-body sector exhibits a quantum anomaly: The RG 
flow of the renormalized coupling Xad exhibits a limit 
cycle, which, due to its periodicity, breaks the classically 
continuous scaling symmetry to the discrete subgroup Z. 
The Efimov parameter can be determined from the pe- 
riod of the limit cycle [23j and is given in our approxi- 
mation by 

,s = « 0.925203. (10) 

The exact result is given by sq ~ 1.00624 Q. Consider- 
ing the simplicity of our pointlikc approximation, which, 
as discussed in Sec. IIII1 does not resolve any momentum 
or frequency structure of the effective (k dependent) in- 
teraction vertex of the three-body sector, the agreement 
is quite good. In fact, in previous work (23j we have 
shown how to obtain the exact value of s using the FRG. 



— |aj Two-body sector 




dt * = dt 





-^b) Three-body sector 




FIG. 2: The flow equations in terms of Feynman diagrams 
for the (a) two-body, (b) three-body, and (c) four-body sec- 
tor. All internal lines denote full, regularized propagators. 
The scale derivative dt on the right hand side of the flow 
equations acts only on the regulators. Solid lines represent 
elementary bosons ip, while dashed lines denote composite 
dimers (f>. The vertices are: Yukawa coupling h (small black 
dot), atom-dimer vertex Xad (open circle), dimer-dimer cou- 
pling X$ (black circle), coupling j3 (two circles), and the atom- 
atom-dimer vertex 7 (black square). Due to the large number 
of diagrams for the latter two vertex functions, we only show 
two exemplary diagrams. 

The presence of N-body bound states leads to diver- 
gencies in the corresponding N-body vertices. The peri- 
odic divergencies in the analytical solution of Xad in Eq. 



- \ c) Four-body sector 
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© correspond therefore to the presence of the infinitely 
many Efimov trimer states at the unitarity point. 

We can use the latter correspondence to calculate the 
bound state spectrum also away from unitarity. The 
trimer binding energies are calculated by determining the 
atom energies for which Xad exhibits divergencies in 
the IR as function of a -1 . The trimer binding energy is 
then given by Et = —ZE^. The result is shown in Fig. [1] 
In this plot the dashed, green lines indicate the binding- 
energies of the Efimov trimer states. For calculational 
purposes we switch to the static trimer approximation 
which is completely equivalent to our two-channel model 
in Eq. @ . We describe this procedure in Appendix [X] 

At the unitarity point the trimer binding energies form 
a geometric spectrum and the ratio between adjacent lev- 
els is given by 



E, 



(n+l) 



E, 



(„) 



(11) 



which can be understood from the limit cycle flow of Xad ■ 
At each scale k = Ae 4 , where Xad diverges, one hits a 
trimer state. The RG scale k can in turn be connected to 
the atom energy 0, [2^, [2!| and as the divergencies 
appear periodically in t one easily obtains Eq. (jlll) . 

There is an additional universal relation obeyed by 
the trimer energy levels which we may take as a mea- 
sure of the quality of our approximation. It is given as 
the relation between the trimer binding energy E* for 
a — > oo and value of a for which the trimer becomes 
degenerate with the atom-dimer (a!j_) and three-atom 
threshold (a*_), respectively. For comparison we define 
a wave number n* by E* = —h 2 K* 2 /M (in our conven- 
tion, E* = -2k* 2 ) and find 



a_K 



-1.68, 



0.08 



[12) 



which has to be compared with the exact result 
a*_n* = -1.56(5), a* + n* = 0.07076 from the fully 
momentum-dependent calculation in 0, [22j ■ The agree- 
ment with our approximate solution suggests that our 
model should provide a solid basis for the step to the 
four-body problem. 



V. THE FOUR-BODY SECTOR 

Recently, the solution of the four-body problem in the 
low-energy limit has gained a lot of interest. In quantum 
mechanical calculations the existence of two tetramer 
(four-body bound) states was conjectured for each of the 
infinitely many Efimov trimers [H,I3|- By calculating the 
lowest few sets of bound state levels von Stecher et al. 
[1, [2(| concluded that both ratios of energies between the 
different tctramcrs and the trimer state approach univer- 
sal constants. However, with the quantum mechanical 
approach the calculation directly at the unitarity point 



(a 



= 0), marked explicitly in Fig. [TJ turns out to 



be difficult, although this point is of great interest when 
one wants to gather evidence for universality of the four- 
body system. In fact, in the three-body sector the infinite 
RG limit cycle appears only exactly at the unitarity point 
and its universal appearance is directly connected to the 
breaking of the continuous scale symmetry. Within our 
approach, the unitarity region is easily accessible. 
In order to investigate the four-body sector we include 
all possible, U(l) symmetric, momentum-independent in- 
teraction couplings in the effective flowing action T^. If 
one assumes all these couplings to be zero at the mi- 
croscopic UV scale A, one can show, by evaluating all 
possible Feynman diag rams and using the vacuum hier- 
archy described in [231 ] . that from all possible four-body 
couplings only the three couplings A , /3, and 7 are built 
up by quantum fluctuations and are therefore included 
in Eq. Couplings other than A^, /3, and 7, such as, 
for instance, the term ~ (i/;*^) 4 are not generated during 
the RG evolution. This consideration leads to our ansatz 
for the effective average action ©. 

For the investigation of the unitarity point we first 
switch to rescaled, dimcnsionless couplings 



A, 



7T 2 /l 4 



fc 4 ^2^,5 

A«>, f3 = —/3, 7=—— 7, (13) 



and obtain the corresponding flow equations by insert- 
ing the effective flowing action Tk, Eq. @, into the 
Wetterich equation ([1]). By the use of the rescaled cou- 
plings we find three coupled ordinary differential equa- 
tions, which are again coupled to the two- and three- 
body sectors, but become explicitly independent of h and 
k. We show the diagrammatic representation of the flow 
equations in Fig. [2jc). Their analytical form at the uni- 
tarity point is given by 3 

128 
125 

111 1 98 

d t X^ = — + -P - -X AD + 3A + -^A 2 (15) 
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3 For illustrative purpose we show the analytical form of the flow 
equations at the unitarity point only. Away from this limit their 
explicit expressions become much more complex. 
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FIG. 3: (Color online) Renormalization group limit cycle behavior of the three- and four-body sector at the unitarity point 
E$ — a -1 = 0. The real parts of the rescaled, dimensionless couplings Xad,i, 4A,£,i, Pi/6, and 71/IOOO are plotted as functions 
of t — ln(fc/A). Not only the three-body coupling Xad,i (dashed, black) exhibits a limit cycle behavior, but also the four-body 
sector couplings A^i (red, dotted), /3i (blue, solid), and 71 (green, solid) obey a limit cycle attached to the three-body sector 
with the same period. 



We pointed out in the last section that the appear- 
ance of bound states is connected with divergent vertex 

(n) 

functions T, and we exploit this behavior to determine 
the bound state spectrum of the three- and four-boson 
system. At this point we must note that these infini- 
ties are complicated to handle in a numerical solution of 
the theory. In particular, the numerical treatment of un- 
bounded limit cycles is problematic due to the periodic 
infinities during the RG flow. In order to circumvent this 
difficulty we used the method of complex extension, de- 
veloped in [HI. The basic idea is to extend the domain 
of the running couplings to the complex plane 

^AD ~ > ^AD.l + i^AD,2 -V — ^ A^l + i\<f, t 2 

On the one hand this effectively doubles the number of 
real flow equations and additional initial conditions must 
to provided. We choose Xad,2 = e = 10~ n in our numer- 
ical calculation and take all other imaginary parts to be 
zero in the UV. On the other hand this procedure allows 
us to perform the numerical integration of the flow equa- 
tions as it regularizes the periodic infinities in the flow 
and makes the numerical treatment feasible. Physically, 
by the complex extension we convert the stable bound 
states into metastable resonances and by taking different 
values of e we are able to vary the decay width of the 
resonances. One may compare this with the procedure 
of Braaten and Hammer [35| who introduce a parameter 
77* in order to model the decay of the trimers to deeply 
bound states which have not been included in the effec- 
tive model. In this line we also view our complex exten- 
sion as a way to include these deeply bound states in the 
FRG calculation. Specifically, we find that for e <C 1 the 
decay width of the n th Efimov trimer rip is given by 

(n) (n) 

r^, = AeEj, at unitarity. This is in agreement with the 
result in Q 

4") « ^4 n) (19) 



which holds for small 77*. Thus, for e <C 1, the relation to 
the parameter 77* introduced by Braaten and Hammer is 
given by 

e=f. (20) 

The result of the numerical calculation of the four-body 
sector at unitarity is shown in Fig. [3] Here, we display 
the RG flows of the real parts of all nonzero three- and 
four-body sector couplings as a function of the RG scale 
t = ln(fc/A). The three-body coupling \ad,i (black 
dotted line) exhibits the well-known limit cycle behavior, 
described in Section ITVl with the period being connected 
to the Efimov parameter sq. Remarkably, there is an 
additional limit cycle in the flow of the four-body sector 
couplings with a periodic structure of exactly the same 
frequency as the three-body sector. This four-body 
sector limit cycle exhibits resonances which are shifted 
with respect to the ones of the three-body system. The 
magnitude of this shift is given by a new universal 
number, which is inherent to the four-body sector. 

Our observation is that the four-body sector is 
intimately connected with the three-body sector at 
the unitarity point. It is permanently attached to 
the running of the three-body sector from the first 
three-body resonance on. From here on the periodic 
structure of the flow remains unchanged as one goes to 
smaller values of k. Due to this tight bond between the 
three- and four-body sector, there stays no room for an 
additional four-body parameter. 

We also find that the magnitude of the shift beyond 
the first resonance is neither dependent on the initial 
values of the four-body sector couplings in the UV nor 
is it influenced by finite range corrections which we 
are able to check by choosing different values for the 
Yukawa coupling h. Arbitrary choices lead to the same 
behavior. Having done this calculation directly at the 
unitarity point our conclusion is, that, within our simple 
approximation, the four-body sector behaves truly 
universal and independent of any four-body parameter 
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FIG. 4: (Color online) Calculation of the universal ratios for 
the lowest five set of levels. The calculation is done for dif- 
ferent values of the Yukawa coupling h 2 /A which determines 
the effective range in our model. The dotted lines are only 
guide for the eye. (a) Ratios between the values of scattering 



lengths CL^lt an d a TVi f° r which the tetramer and correspond- 
ing trimer become degenerate with the four-atom threshold, 
(b) Ratios between the values of binding energies E^ t and 
^Tri a ^ resonance a — > oo. 



confirming the conjecture made by Platter et al. and 
von Stecher et al.. We expect that universality will also 
hold for an improved truncation. 

Naively one expects that each resonance in the flow 
of the vertex functions is connected to the presence of a 
bound state. As one observes there are also additional 
resonances in the four-body sector being degenerate with 
the three-body sector resonances. However, we arrive at 
the conclusion that these resonances are artifacts of our 
approximation. The mathematical structure of the flow 
equations is of a kind that divergencies in the three-body 
sector directly lead to a divergent four-body sector. We 
are confident that the resonances at these positions will 
disappear as one includes further momentum dependen- 
cies in the field theoretical model. Therefore we can al- 
ready infer from the calculation at unitarity that within 
our approximation we are only able to resolve a single 
tetramer state attached to each trimer state also away 



from unitarity. In contrast, the "exact" quantum me- 
chanical calculations in [l)-^, [2(| predict the existence of 
two tetramer states which have recently been observed 
by Ferlaino et al. [2l[. As one includes further momen- 
tum dependencies, it is well possible that not only the 
degenerate resonance disappears but also new, genuine 
resonances associated with the "missing" tetramer state 
will appear at the same time. This effect indeed occurs 
in the three-body problem. There, it is essential to in- 
clude the momentum dependent two-atom vertex. Only 
under this condition one arrives at the quadratic equa- 
tion as in © which gives rise to the Efimov effect. This 
can easily be seen by taking a look at the flow equation 
of Xad depicted as Fcynman diagrams in Fig. &b). The 
assumption of a momentum independent two-atom in- 
teraction corresponds to a momentum (and frequency) 
independent dimer propagator. In this approximation 
the first term on the RHS of Fig. Hfb) vanishes because 
all poles of the loop frequency integration lie on the same 
complex frequency half-plane. This directly leads to the 
loss of the Efimov effect in this crude level of approxima- 
tion. 

We can also use our model to investigate the full bound 
state energy spectrum by solving the flow equations for 
arbitrary values of the scattering length a. The energy 
levels of the various bound states are then determined 
by varying the energy of the fundamental atoms such 
that one finds a resonant four-body coupling in the IR. 
The result of this calculation, using the static trimer ap- 
proximation presented in Appendix is shown in Fig. 
[TJ where we plot the energy levels of the various bound 
states versus the inverse scattering length. We find one 
tetramer state attached to each of the Efimov trimer 
states. These tetramer states become degenerate with 
the four-atom threshold for negative scattering length 
and merge into the dimer-dimer threshold for positive 
a. In the experiment this leads to the measured reso- 
nance peaks in the four-body loss coefficient. In order 
not to overload the plot we show only the first three sets 
of levels, although the FRG method allows to calculate 
an arbitrary number of them. One also observes that 
the shape of the tetramer levels follows the shape of the 
trimer levels. In analogy to the three-body sector one can 
calculate a universal formula relating a tetramer binding 
energy E^ et — — 2k^ 2 at a — > oo with the corresponding 
scattering length at which the tetramer becomes degen- 
erate with the four-atom threshold and the dimer- 
dimer threshold a^ + , respectively. We find 



-1.75, 



0.20. 



(21) 



In their recent quantum mechanical calculations von 
Stecher et al. were able to calculate the lowest few sets of 
bound state energy levels [|| . From their behavior it was 
inferred that the ratio between the tetramer and trimer 
binding energies approaches a universal number within 
these first few sets of levels. Figuratively speaking it is 
therefore expected that the universal regime in the en- 
ergy plot in Fig. [1] is reached very fast as one goes to 
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smaller a 1 and E^. 

In order to investigate this observation we calculate 
the behavior of two ratios as a function of the set of level 
for which they are determined. The first ratio relates 
the negative scattering lengths a^ t and a^ft f° r which 
the n th tetramer and trimer become degenerate with the 
four-atom threshold. The second is the ratio between the 
binding energies of the n th tetramer E^l and the n th 
(n) 

trimer Ej, ri at resonance, a — > oo. The resulting plots 
are shown in Fig. 21 Wc calculate the ratios for different 
values of the microscopic couplings in order to test the 
degree of universality of the various sets of energy levels. 
In the plots we show in particular the dependence on the 
choice of the Yukawa coupling h determining the effective 
range r g of the model. As one sees, only the first of the 
ratios depend on the microscopic details. Already from 
the second set of levels on the microscopic details are 
washed out and the ratios become independent of the 
choice of initial conditions: The regime of universality is 
reached extremely fast and as a -1 and E^ are lowered one 
will ultimately hnd the four-body limit cycle described 
above. 

For the asymptotic ratios we find 

a™ « 0.518 4ft, (22) 
4"j ps 4.0174rt- (23) 

Von Stecher et al. find a^/a^ft ps 0.43(0.9) 
for the deeper (shallower) bound tetramer and 

E^ljE^ ps 4.58(1.01), respectively. Considering 
the simplicity of our model the agreement is quite good. 
With an ultracold bosonic Cs gas Ferlaino et al. found 
a^l/a^ri ~ 0.47(0.84). In this remarkable experiment 
only the lowest set of tetramer states in the energy 
spectrum had been accessible due to the particular 
scattering length profile. Considering our observation 
that the deepest set of levels is still strongly dependent 
on the microscopical details it cannot be expected to 
find the universal numbers in this particular setting. 
Therefore more experiments for bosons interacting via a 
larger scattering lengths would be desirable. 



VI. CONCLUSIONS 

In this paper, we investigated the four-body problem 
with the help of the functional rcnormalization group. 
Employing a simple two-channel model with pointlikc 
three- and four-body interactions we were able to inves- 
tigate universal properties at the unitarity point a — > oo, 
E^, = as well as to perform computations away from it. 
In the RG language the Efimov physics of the three-body 
problem manifests itself as an infinite RG limit cycle be- 
havior of the three-body coupling constant at unitarity. 
We found that also the four-body sector is governed by 



such a limit cycle which is solely induced by the RG run- 
ning of the three-body sector, signaling the absence of a 
four-body parameter. 

We also computed the energy spectrum away from uni- 
tarity and were able to obtain the universal relations be- 
tween four- and three-body observables in our approxi- 
mation. Our calculation provides an explanation for the 
findings of von Stecher et al. who found that these 
ratios approach universal constants very quickly as they 
are computed for higher and higher excited states. We 
also found a dependence of the ratios for the lowest level 
on microscopic details such as the effective range. This in 
turn is of relevance for the experimental observations by 
Ferlaino et al. 21]. In this experiment the lowest states 
have been measured and one can therefore not expect to 
find the exact universal relations between them. 
Considering the simplicity of our model, the agreement 
with the previous studies in [l|-[3| is quite good. There 
had been some disagreement in literature about univer- 
sality and the absence or existence of a four-body pa- 
rameter, see e.g. (36l - [38| . Our RG results support the 
conclusion that the four-body system is universal and in- 
dependent of any four-body parameter. 
An important shortcoming of the pointlike approxima- 
tion is the absence of the shallower of the two tetramer 
states. Obviously the pointlike approximation of the 
three- and four-body sectors is not sufficient and in fu- 
ture work one should include momentum dependent in- 
teractions. From the energy spectrum in Fig. [T] it be- 
comes also evident that the excited tetramer states can 
decay into an energetically lower lying trimer plus atom. 
The higher excited states in the four-body system are 
therefore expected to have an intrinsic finite decay width 
0. Whether this width has a universal character still 
remains an opened question as well as in which way the 
corresponding imaginary coupling constants will change 
the RG analysis. 

The inclusion of the full momentum dependencies in the 
three- and four-body sector seems to be a rather compli- 
cated task. In the effective field theory study of the three- 
boson system the introduction of a dynamical dimer field, 
often called the di-atom trick Q , has been a decisive step 
towards the exact solution of the three-body problem. 
From this perspective we suggest that the inclusion of a 
dynamical trimer field in the effective action might help 
to simplify the momentum dependent calculation. 
The four boson system remains still a subject with many 
open questions. With our RG analysis in the point- 
like approximation, we made the first step towards a 
rcnormalization group description of the four-body prob- 
lem supplementing the previous quantum mechanical ap- 
proaches. From this perspective this work provides a 
starting point for a deeper understanding of universality 
in the four-body problem. 
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Appendix A: The trimer approximation and 
rebosonization 



In this appendix we will apply the rebosonization 
method developed in [4(| to our model ©. We intro- 
duce an additional trimer field x, representing the bound 
state of three bosons, which then mediates the atom- 
dimer interaction. A similar procedure had already been 
used long time ago by Fonseca and Shanley [TtJ in the 
context of nuclear physics and was recently employed by 
us in [2^, [29| for the treatment of the three-component 
Fermi gas. There are several reasons for employing this 
procedure. First, it is useful to reduce the number of res- 
onances one has to integrate through in the RG flow. In- 
stead of calculating the divergent coupling Xad one only 
has to calculate zero-crossings of the trimer energy gap 
which is numerically much easier to handle. Secondly, by 
the introduction of a dynamical trimer field one may be 
able to mimic some of the complicated momentum struc- 
ture of the atom-dimer interaction in a simple way which 
could probably be sufficient to find the missing tetramcr 
state in our calculation. The third point is of a more tech- 
nical nature and concerns the method of rebosonization, 
which we will employ here in a quite extensive manner. 

In the three-body sector the real atom-dimer coupling 
Xad exhibits divergencies when the energy gap of the 
fundamental atoms E^, is tuned such that one hits the 
trimer bound state in the IR. In the static trimer approx- 
imation, the coupling Xad is mediated by the exchange 
of a trimer field \ with the non-dynamical, inverse prop- 
agator P x = mt, which can be depicted as 




X 




A, 



+~(<W + + Xad4>*^*H 

+g{x*W> + x<A*V>*) 

+5 3 { x *rw^ + xu*rr)}- (ai) 

The Yukawa interaction g couples the trimer field to the 
dimer and atom field. The Si are the additional U(l) 
symmetric four-body couplings which are generated by 
quantum fluctuations. All other possible couplings can 
be shown to stay zero during the RG evolution provided 
they are zero at the UV scale. Also the coupling Xad 
is regenerated through a box diagram in the RG flow. 
However, it is possible to absorb all these emerging cou- 
plings by the use of the rebosonization procedure. 
For this matter we promote the trimer field x to be ex- 
plicitly scale dependent, x ~ ► Xk, X* ~ * Xk an< ^ the Wet- 
terich equation generalizes to 



d k r k [$ k } = ^Tr (ri 2) [$ fe ] 



d k R k 



5<P, 



■T k [* k ])d k $ 



(A2) 



where $fc now includes all fields including the trimer 
fields (x, X*)- The additional term in the generalized flow 
equation (|A2[) allows for the absorption of the reemerging 
couplings since one has the freedom to choose the scale 
dependence of the trimer fields as a function of fields. In 
order to continuously eliminate the couplings Xad and Si 
wc choose 



dkXk 
dkXk 



H>Ca,k + tp*Xktp(b,k 

+u*cp*Cc,k + u*rrUk- 



(A3) 



Upon inserting Eq. (| A3|) into the generalized Wet- 
terich equation (| A2|) the condition that the flows of Xad 
and Sj vanish leads to 



Cc 



d k X 



k*AD 



2ff 

dkS 2 



Cd = - 



d k 5i 
2ml 

dkS 3 + g(b 



(A4) 



The trimer field x ~ "0 3 is introduced on the microscopic 
scale by a Hubbard-Stratonovich transformation and our 
ansatz for the effective average action, motivated by the 



When one calculates now the flow equations of the re- 
maining flowing couplings by projecting Eq. (|A2[) onto 
them, one obtains new contributions due to the presence 
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of additional terms arising from Eq. (| A3|) . 
In our static trimer approximation the trimer field has no 
dynamical propagator and the model given by Eq. (|Alj) 
is completely equivalent to the two-channel model in Eq. 
([2]). Furthermore, no regulator has to be specified for 
the trimer field since in our approximation the original 
atom-dimer coupling Xad is solely replaced by g 2 /rn 2 ^. 
At this point it already becomes clear why the modified 
flow equations will be easier to handle numerically: in- 
stead of calculating a divergent Xad hi the three-body 
sector one has only to deal with zero crossings of m£ at 
the values of k where originally Xad had divergencies. 
The modified flow equations are given by 

dtg = d t g\® k +m 2 x ( a , 

da = 9*7k + 2sCi, (A5) 



where the first terms in the flows are the original flow 
equations with the trimer field taken to be scale indepen- 
dent. In fact, by expressing all flow equations in terms 
of the coupling G = g 2 one can also get rid off the prob- 
lematic g in the denominator of Ca in Eq. (|A4[) . We 
point out that the static trimer approximation allows to 
calculate easily the three-body sector. In the four-body 
sector the original divergencies of Xad still appear since 
trimers \ appear in the corresponding flow diagrams and 
therefore one has to deal with terms ~ V m x' ^ or * n * s 
reason it had been essential to perform the complex ex- 
tension described in Section [V] Finally, the bound state 
spectrum of the trimers can be computed by calculating 
the poles of the trimer propagator, m^k = 0,E^) = 0, 
in a straightforward manner. 
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